Form factors of twist fields in the lattice Dirac 
theory 



o 

(N 



dt 



X 



P Gavrylenko 1 , N Iorgov 2 and O Lisovyy 3 

1 Department of Physics, Kyiv National University, 03022 Kyiv, Ukraine 

2 Bogolyubov Institute for Theoretical Physics, 03680 Kyiv, Ukraine 



^ ■ Laboratoire de Mathematiques et Physique Theorique CNRS/UMR 6083, 



'. Universite de Tours, Pare de Grandmont, 37200 Tours, France 

\Q . E-mail: iorgovObitp .kiev. ua, lisovyiOlmpt .univ-tours . f r 

Abstract. We study U(l) twist fields in a two-dimensional lattice theory of 
massive Dirac fermions. Factorized formulas for finite-lattice form factors of these 
fields are derived using elliptic paramctrization of the spectral curve of the model, 
iQ} ' elliptic determinant identities and theta functional interpolation. We also investigate 

i-C . the thermodynamic and the infinite- volume scaling limit, where the corresponding 

expressions reduce to form factors of the exponential fields of the sine-Gordon model 
at the frcc-fermion point. 
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1. Introduction 

It is a general property of two-dimensional quantum field theories that their symmetries 
give rise to new local fields, whose correlation functions are nontrivial even if the 
underlying theory is free. Two paradigmatic examples are given by the disorder variables 
in the Ising field theory [16] 130] and U(l) twist fields in the massive Dirac theory 
[23| [29| |3T] , directly related to the exponential fields in the sine-Gordon model at the 
free-fermion point. 

Correlation functions of twist fields in the Dirac theory, as well as in its 
generalizations to curved space and non-zero background magnetic field [6j [7J HEJ EZ] , 
are also interesting from the mathematical point of view. They satisfy nonlinear 
differential equations [U [HI |29] , which in the simplest cases can be solved in terms of 
Painleve functions. The knowledge of the long- and short-distance behaviour of the two- 
point correlators provides solutions to nontrivial asymptotic and connection problems of 
Painleve theory. Recently, it has also been observed pjJ] that such correlators coincide 
with the gap probabilities for the classical kernels arising in the representation theory 
of big groups [2] . 

The aim of this paper is to construct lattice analogs of U(l) twist fields in the Dirac 
model, satisfying the following properties: (i) they should be defined via the branching 
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of lattice fermion fields, (ii) one should be able to calculate their form factors explicitly 
and (iii) these form factors should reproduce the known expressions in the scaling limit. 
Besides full control of the theory, such an integrable finite-lattice regularization can 
be used for investigations at non-zero temperature and for a mathematically sound 
derivation of the relative normalization of conformal and infrared asymptotics of the 
two-point correlator [51 [21]. It may also be instrumental in going beyond the free- 
fermion point. 

While the Ising field theory possesses a natural lattice regularization, only a few 
results are available in the Dirac case. First attempt to introduce twist fields on 
the infinite lattice was made in [25]. The corresponding definition was supported by 
the computation of the vacuum expectation value (reproducing the expected scaling 
dimension), and was further strengthened by the analysis of correlations at the critical 
point [26]. Another, seemingly unrelated definition was used in j5j [20] to derive a 
number of determinant representations for the two-point function of lattice twist fields. 
The present work is devoted to the computation of their form factors, i. e. matrix 
elements of the field operators in the basis of transfer matrix eigenstates. 

The paper is planned as follows. In Subsections 12.11 and 12.21 we introduce a one- 
parameter generalization of the lattice Dirac operator considered in [51 [20] and explain 
the definition of twist fields in the functional integral framework. Transition to the 
operator formalism is performed in the next subsection. Using the coherent states 
approach, the transfer matrix and the twist field operator are written as exponentials 
of fermion bilinears, see formulas (jHJ) and ( ITT]) below. Subsection 12.41 is devoted to 
the construction of multiparticle Fock states simultaneously diagonalizing the transfer 
matrix and the operator of translations. 

In Subsection 13.11 it is explained how form factors of twist fields can be found from 
the linear transformations relating fermionic creation-annihilation operators of different 
periodicity. In particular, the vacuum expectation value and two-particle form factors 
are expressed in terms of two square matrices C and D of dimension equal to the 
lattice size. Essentially, one needs to compute the quantities D~ x , D~ X C and det.D, 
cf. (j2"S"l) - (j3"2"|) . This task is solved in Subsection 13.21 by first noting that in the elliptic 
parametrization of the spectral curve of the model C and D are given (up to diagonal 
factors) by elliptic Cauchy matrices, and then using Frobenius determinant identity and 
theta functional interpolation along the lines of [13]. The corresponding expressions are 
further simplified in Subsection 13.31 Finite-lattice two-particle form factors are given 
by (|42p . (jUJ) and (I45p . and the multiparticle ones have the factorized form fT46l) . fj49|) . 
These formulas represent the main result of the paper. In Subsection 13.41 we analyze 
the thermodynamic (infinite- lattice) limit. The final answer has a remarkably simple 
expression in terms of the Jacobi theta functions, see (I54|) - (l56|) . We remark that the 
vacuum expectation value (}54"1) reproduces the earlier result of [25]. Field theory limit 
is considered in Subsection 13.51 It is shown that the scaled form factors coincide with 
those of the exponential fields of sine-Gordon model at the free-fermion point [U [231 El] • 
We conclude with a brief discussion of results and open problems. 
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2. Lattice Dirac theory 

2.1. Fermions 

Let if), if) denote two 2-component Grassmann fields onanMxJV square lattice. Consider 
the standard free-fermion action S[il>, ip] = if)Dif), where the lattice Dirac operator is 
chosen as 




s x c y^y c y ~l~ c x^x 

c x^ '—x c x s y ~ s x c y^ ' -y 



(1) 





Here V x ,y denote the shifts by one lattice site in the horizontal and vertical directions, 
so that e. g. V x if) XtV = if) x+ i jy , V y if) XtV = ift x ,y+i- The boundary conditions with respect 
to x and y are antiperiodic and a-periodic, respectively. This means that 

~ ^ **Px,y 

The parameters q, Sj, c*, s* {i = x, y) are expressed in terms of two constants JC XtV G M>o 
as 

Q = cosh2/Q, Sj = sinh2/Cj, c* = cosh 2/C*, s* = sinh 2/C*, 

where the dual couplings are given by )C* = arctanh e~ 2ICi . They satisfy Ising-type 
relations s* = s" 1 , c* = CjS^ 1 . It will be assumed in the following that /C* < JC y . Dirac 
operator considered in [20] is obtained from (CQ) by setting K x = fc y . 

Fermion propagator can be found using Fourier transform. One obtains 

1 p iH x - x ')+ ie (y-y') ( r * o — q* r p~ ie r — r*p^ \ 

Wx,vH>*,vl MN Z-^2c y (cosh 7fl - cos 4>) \ -c y + c* x e-^ c* x s y - s x c y e lf) ) K) 
where 

cosh jg = c* x c y — s* x s y cos 6. (3) 

The multipoint correlations are readily computable from the Wick theorem. The 
summation in f^j is performed over = || (j + |), j — 0, . . . , M — 1 and 9 = ^ (k + a), 

k = 0, . . . , N — 1. In the thermodynamic limit M,N — > oo one has jt^^2 — > 

<p,e 

i r r 27r 



ff Q n d(j) d6 and the answer can be expressed in terms of elliptic integrals. 



2.2. Twist fields 

It is instructive to start with an example. Choose a closed path V on the dual lattice 
(Fig. 1) and make the transformation if) i— > e 2nlu if), if) i— > e~ 27TlI/ if) with i/Glat all lattice 
sites inside this contour. Because of the global £7(l)-symmetry, the action S[if), if)] will 
change only at the edges intersected by V. 

If the corresponding changes are made along an open path Vab joining two points A 
and B on the dual lattice, the resulting functional integral will depend on the positions 
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of these points and homotopy class of the path, but not on its shape. One way to choose 
Vab is shown in Fig. 2. In this case, the action is modified by 



SSab = SSa + SSb + SS\ 



b.c. 



(4) 



where 5Sa,b correspond to the vertical segments and <5»Sb.c. to the horizontal one. 
Explicitly, 

ii-i 

SS A = 2* sin^ ( eiWU rt-i,y»<y» + ^Rtf'ti-w) . 

•y"=0 

(J5 fl = -2zsin™ £ (e^^^^e,/ + , 

y"=0 



x'-l 



5S h . c . = -2ismnv^- ^ (e^+^V^Ar-iV^o - 



«r(a+a') 72 -2 



i J x",oi ) x>'N-l 



with a' = a + za Without any loss of generality, we assume that < a, a' < 1. 

Twist fields live on the dual lattice. Their two-point correlation function is defined 
as the normalized partition function of the Dirac theory with the defect contribution 



(O a>a ,(A)O a , >a (B)) 



S&M+SSab 



(5) 



/ VifiVtp e s ^ 

Note that the effect of 5S^ C _ amounts to the change of the vertical boundary conditions 
for fermions from the horizontal interval [x, x' — 1] from a- to a'-periodic ones. More 
generally, the appearance of twist field O a , a '{A) in an arbitrary correlation function 
means that 

• the vertical boundary conditions for fermions to the left and right of A are a- and 
a'-periodic, respectively; 
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• the term 5Sa should be added to the action. 

In the heuristic continuum limit, this corresponds to integrating over field configurations 
having counterclockwise monodromy e 2mv of if) (resp. e~ 2mu for ip) around A. 



2. 3. Coherent states and the operator formalism 

Locality of twist fields becomes manifest in the operator formalism, which also provides 
a convenient framework for the computation of correlation functions. 

Let us introduce two sets of fermionic creation-annihilation operators satisfying 
canonical anticommutation relations 

{a y , a^} = {b y , b ] yl } = S yy >, y, y' = 0, . . . , N - 1, 

with all other anticommutators vanishing. Define in the usual way the vacuum vectors 
(vac\ and \vac), normalized as (vac\vac) = 1, and the corresponding 2 2Ar -dimensional 
Fock space JF '. Further, introduce the coherent states 

N-l 

= ex P Yl ( a l^iy + 6 W,y)l«ac), 



y=o 



N-l 



V>J = (vac\ exp i^l,y a V + ^,A)> 

y=0 

where {i^ l x>y }, {i^x^} denote Grassmann variables anticommuting with all creation- 
annihilation operators. These states satisfy the following standard properties: 

• For y = 0, . . . , N — 1 one has 

a y |^i,^>=^ v |^i,^>, 
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N-1 

M,^> 2 Ml, j>l) = eX P (-^iv^x'.y + ^l'^l,y) ■ ( 6 ) 

y=o 

The identity operator can be represented as a 4iV-fold Grassmann integral 



J v^lv^lv^v^l faMtfA x (7) 

N-1 

x exp $lv1&, v ~ $Ul,y) ■ 

y=0 

The trace of any operator O can be written as an integral of its matrix element in 
the basis of coherent states with a Gaussian kernel: 

Tr O = J V^V^V^V^ MrfllOWl, x (8) 

N-1 



N-1 
y=0 

.1 ."7. 1 ' 



Matrix elements (ij) x ,,ip x ,\0\xj) x ,ip x ) can be obtained by writing O in normally 
ordered form, making therein the replacements 

and multiplying the result by (jBJ). 
Now consider the operator 



N-1 



V a = : exp - c x x c y V-y)a\ + (c* ^ - 1) (aja y + + 

y=o 

+ a y( s y- c x lc y^y) b y} : > ( 9 ) 



where it is understood that 

l y+N ) = e 2nia ( a \ 

°y+ N J \ by J \ U y +N J \ u y 

Rewrite the quantity Z = Tr V^f using at the first stage (ED with x = N — 1, x' = to 
calculate the trace, and inserting M — 1 resolutions of unity (J7J) (with x = A; — 1 and 
x' = k between the k-th and (k + l)-th factor of V a ). Then, computing matrix elements 
of V a in the basis of coherent states, the reader may easily check that Z coincides with 
the partition function of the Dirac theory described by ([!]). 

In fact V a is the transfer matrix characterizing discrete time evolution of twist fields. 
Vertical defects of the action associated to them divide the horizontal axis into intervals. 
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The evolution in different intervals is governed by the matrices V a with appropriate 
values of a. Twist fields are represented by the operators 

2/-1 

O a , a > (y*) = exp lixiv ^ (- a l' a y' + b l> V) > ( n ) 

y'=0 

where y* = y — \- This can be seen by noticing that 

N-i „ 
exp{6S A + £ ('•; i,/' --/ - I,/)} = / (e a ,?|0«^ (y*) I^M 1 ) x 

N-1 

X exp {$x-l,y>£yi ~ ^y'^x-l^' + Cy'Va.y' ~ ^x^'^y') > 

and repeating the procedure used above for the computation of Z. For example, two- 
point correlator (jSJ) can be expressed as 

(O a y{A)O a/>a {B)) = Z-'Ti (o a>a , (y*) V$~ x O a ,, a (y'*) i/f . 

The problem of effective calculation of correlation functions of twist fields therefore 
reduces to the computation of form factors of the operator ([IT]) between the eigenstates 
of V a and V a >. Observe that e. g. O aja > (0*) is given by the identity operator on T . 
However, its form factors are nontrivial since the corresponding bra and ket states 
diagonalize different transfer matrices. 

2.4. Transfer matrix diagonalization 

Define Fourier transforms of the creation-annihilation operators: 



t \ ! JV-i / + 
ai 1 1 / a 



-T, 



y 1 e ~ ie y 



!J= 

N-1 



a v \ e idy 



where 6 belongs to the set 6 a = {jr{k + a) \ k — 0, . . . , N — l}. The only nonvanishing 
anticommutators are given by {a^aj,} = {bg,bg,} = dee 1 - The transfer matrix V a is 
block-diagonal in the Fourier basis, 

V a = : exp { ( s y - c x 1 c y e~ %0 ) b\a\ + [c*~ x c y - 1) (a\a e + b\b^j + 
eee a 

+ { s y ~ c x lc y el6 ) a eh 
The conjugation of fermions by V a induces linear transformations 

v„ ( f t J v- 1 = MS) 1 al 
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where A{6) is a Hermitian matrix with unit determinant, explicitly given by 



A(6) 



c y V c x s y s x c y e c x 



U<(9), by a 

e 16 I 

unitary transformation. Here 70 is defined as the positive solution of ([2]), and the 

- ] are the eigenvectors of A(6>) normalized so that 
h{&) h{&) J 

\fi(@)\ 2 + \h{@)\ 2 = 1- The freedom in the choice of the phase can be used to set 

h{d) \ 1 / e K 'y 
h{9) j 2^ 1 e- K 



"y 




(12) 



where 



Xe = [X-e\ 



(1 -ae ie )il -(3e~ ie ) 



(13) 



1 -/3e ie )(l -ae- ie ) 

and a = tanh/C* coth/C^, /3 = tanh/C* tanh/C y . Note that under the above conventions 
one has (3 < a < 1. Root function in (fTB"]) is taken on the principal branch. 
A new set of the creation-annihilation operators 

satisfies canonical anticommutation relations {c^cj,} = {do,d\,} = 5ggi (the other 
anticommutators being equal to zero) and diagonalizes the transfer matrix V a : 

N 



° x ' eeo a 



Introduce the vacua a (yac\ and \vac) a annihilated by all {c e }, {d e } and {eg}, {dg}, 
respectively, and normalized as a (vac\vac) a = 1. Left and right eigenvectors of V a are 
then given by the multiparticle Fock states 

0,(0; (f>\ = a (vac\c0 1 ...c 0m d < j, 1 ...d < j >n , (14) 
\0;(f)) a = c\ i ...c\ m d\ l ...d\ n \vac) a , (15) 



and the corresponding eigenvalue is equal to ( ^ J exp< Id ~ 2~2 le\- 

The states (jl4p - (H5p simultaneously diagonalize the U (l)-charge and the translation 
operator, 

N-l 

Q=Y1 { c Q Cd ~ d e dd ) = { a d ad ~ b e be ) = ( a l a y ~ b l b y) ' 
eee a eeo a y=o 

N-l 

T a = exp ^2 W (-<4 c e + 4°^) = : exp ^ (a y a y+1 - a y a y + b y b y+ i - b\b y ) : . 
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The latter satisfies, e. g., 




under boundary conditions f fTUj) . 

It is useful to note that the twist field (fTTj) is ?7(l)-neutral and coincides with the 
identity operator twisted by translations of different periodicity: 

O a , a , {y*) = [0*>, a {y*T l = T-yT y al . 

This can be seen by comparing linear transformations of fermions induced by the left and 
right side of this relation, and using that their action leaves the vector \vac) invariant. 
Form factors of O a ^ (y*) are therefore given, up to simple multiplicative factors, by the 
scalar products of the Fock states (|T4|) - (|T5|) : 

a (0;#|cV (?/*)|0';</>V =ex P u/( °~ Yl B ) *(0;<l>\o';<f>')a>, 

eeeucp' ee4>ue' 

where 0,0 C 9 a and 0' , <f> C Q a >. Such scalar products normalized by the product of 
the vacua will be denoted by 



T a>a ,{0-<f>\0'-^) = 



a (0;0|0 Q , Q ' (0*) 


0' ; (f)') a > 


a (vac 


O a ,a' (0*) 


vac) a i 



(16) 



3. Form factors 



3.1. General setting 

Let us first recall a few results from (TTJ [131 EE]- Consider two sets of 2L fermionic 
creation-annihilation operators generating equivalent Fock representations in the same 
space. Denote the corresponding 2 L x 2 L matrices by {^i}, {Vi} and {^}, {<p\} with 
i = 1, . . . , L and combine them into L-columns i/j, ip\ cp, cp\ Suppose there exists a 
unitary operator a such that 

V? f A _j / A B \ f ^ 



(p )° I C D IU 

where A, B, C, D are some L x L matrices. The unitarity of a and canonical 
anticommutation relations imply that B = C, A = D and 

DC T + CD T = 0, DB-t + CC^l. 

Suppose that D is invertible. Then, up to inessential phase factor related to the 
choice of the vacua, one has 

^,(vac\o-\vac) v = |det D| 2 . (17) 

General matrix elements of a between Fock states of different types can be expressed as 

^(vaclifji, . . .^i m o^ h . . . (pijjvac)^ = |detD|5 . pf (18) 
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where the entries of the blocks of (m + n) x (m + n) skew-symmetric matrix R are given 
by the normalized two-particle form factors 

(R IXI ) kl = (D^C) . fc . , M = l,...,m, (20) 

(i2jxj) H = -(«Jx/) Ifc = D i ;;., fc = l,...,m, Z = l,...,n, (21) 
(J2 JXJ ) H = (CD- 1 ) . u . , fc,i = l,...,n. (22) 

In the case of interest here, a is the identity operator and L = 2N. The 
creation-annihilation operators in each set are labeled by their £/(l)-charges and the 

/ c 

corresponding momenta. Thus e. g. ib and <p are given by the 2iV-columns 

v d 

built from the operators eg, dg with 9 G a and 9 G a >, respectively. A, B, C, D can 
therefore be seen as block 2x2 matrices with block entries indexed by 9 G a >, 9' G 6 a . 
To find their explicit form, note that for 9 G 6 a i one has 

~t \ 1 . 1 _ p -2-Kiv I t 

W) I - 1 \ r X C 7-T+ / n\ TT/nl\ I 



]v E T^iy («'> X ■ (23) 

* E wo ( | ) • <*) 

where is defined by (fl2l) . Therefore, introducing the notation 

V = e ! ^ 2 5,,, 9,9' eO a , 

K,e' = e l{nu - 9)/2 5 e , e ,, 9,9' eO a ,, 

_smnv f 1 (9)f 2 (9')-f 2 (9)f 1 (9') 

^6,0' — 77 i 1)!— a , " G O a >,V G (7 Q , 

iv sin 

U Qfi' — 77 ; ai—a ) f t fa'jf t U a, 

iV sm 2 -jp 

we find that 

A'CA \ _ / -A'DA 

-A'CA J ' ~ I -A'DA 

This in turn implies that the vacuum expectation value of twist field and its non- 
zero two-particle form factors are given by 

a (vac\O ata '(0*)\vac) a i = |det-D|, (25) 

F a , a >(9;\9';)= -e^-^D^,, (26) 

JV(;0|;0')= -e^'+^D-y (27) 

Fa<*{0\ff\\)= -e^'^ 2 {D~ l C) eei , (28) 

F a , a '(;\e;9')= -e^' 2 {CD' l ) eei . (29) 
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One also has 



Fa,a'{\ Q') = —Fa',ot{9] 9'\] ), (30) 

although this is not immediately obvious from (I28l)-(l29l). This reduces our task to the 
computation of determinant and inverse of D and of the product D~ X C . 

The relations (l23"]) - (!24p imply that the elements of C and D remain invariant if one 
multiplies U(9) from the left by a unitary matrix independent of 9. Together with (1T2"j) . 
this gives 

sin tip XeX-e' ~ X-eXe> 
ziv sm 

sin 7rv XeX-e' + , QoA 
ziv sm 

The matrices C and D are therefore very simply related to the ones appearing in the 
Ising model theory, see Lemma 3.3 in [13]; the main and almost only difference between 
the two cases is the change in the spectrum of quasimomenta. We will now follow [T2j 
to establish elliptic representations for C and D, which then will be used to calculate 
det£>, IT 1 and D~ l C. 

3.2. Elliptic parametrization 

Spectral curve ([3]) is a torus which can be conveniently uniformized by the Jacobi elliptic 

sinh 2K* 

functions of modulus k = — - — —r-. Let us denote K = K(fc), K' = K(\/l — k 2 ), where 

sinh 2/Cy 

K(&) stands for the complete elliptic integral of the first kind. Then [13j[2l] the functions 

sn(u + in) 

z ( u ) = —7 r-r, (33) 

sn(u — ir]) 

X(u) = [ksn(u + irj) sn(u — irj)} -1 , (34) 

with r\ G (— ^-,0) determined by sinh2/C x = isn2ir] } satisfy the relation ([3]) written in 
the form 

s x (A + A^ 1 ) + s y (z + z^ 1 ) = 2c x c y . 

The formulas ( J33i) -(l34 j) bijectively map the real interval C u = {u\R,eu G 

[-K,K),lmu = 0} to C e = {(z, A) = (e ie ,e^)\ 9 G [0,2tt)}. The inverse image of the 

point (e l6, ,e 7fl ) G Cq in C u will be denoted by u$. Note that for 9 G (0,7r] one has 

uq = —u 27T ~0. It is also useful to define the function xq = — 7, which continuously 

2K 

increases from — | to | when 9 varies from to 2n. 

Lemma 4.1 in [13] shows that under the above parametrization matrix elements 
(|3T]) -( l32|) can be written as 

— is* sin7rz/ 

Ce,e> = AT , . f .1 cn ( u e + u e >), (35) 
A* v sinh 70 sinh 70/ 

-SySinvrz/ dn(u e - u g >) , . 

Uq 0i = = . [oO ) 

sinh 70 sinh 70/ sn(ug — uq') 
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We will also need Jacobi theta functions $i..a{z) of nome q = e l 
to the elliptic modulus and half-periods by 



which are related 



IK 



UK' 



where $i = $i(0) for i = 2, 3, 4. The quantities det D and D" 1 can be computed in terms 
of these functions using Frobenius determinant identity, while D~ 1 C can be found from 
a theta functional analog of the Lagrange interpolation formula (see Section 5 of [13] 
for the details of a similar cumbersome calculation). The result is as follows: 



detD 



MX a 



s y sin rev $2$4: 



N 



? 



N 



-1 



A' 



J sinh 70 ] [ sinh 



x 



x ■ 



eee a eao a , 

Yle,e>eo a ,e<e> M x e> ~ x e) Ue ,e'eo„,,e<fl' M x e - x B >) 



Uee0 a ,e>eo al M x 0' ~ x o) 

M X a' - X a + X 6 -Xff) 



(37) 



s y sin nu "&2$i 



[sinh 7 9 sinh jgi 



IL» e o„, M x e ~ x &") Yle»eo a M x &' ~ x e") 

x B ») 



2 x 



I\e"ee a ,6"^eM x 6 - x e >>) Yle"ee a ,,e"^e' M x e' 



(38) 



—i- 



. #2 ( X a > - X a + x e + xg/) 



sinh 70 ] I 
sinh 70/ 



x 



Yl 9 »e9 a , M x e - x 6 ") n.o»eo a ,o»#) M x e> 



x d ". 



(39) 

Ue"ee a/ M x 0' + x 6") IU»e8 a ,e»& M x o ~ x e >>) ' 
with X a = J2eee a X Q- ^ n ^ ne nex ^ subsection, the answer (p7jl - (j3*9j) is rewritten in a 
somewhat different form, which turns out to be more suitable for the analysis of the 
thermodynamic limit and for the computation of multiparticle form factors. 



3. 3. VEV and form factors 

We illustrate the procedure by transforming the expression for (D~ 1 C) e e ,. First rewrite 
the second line of (1391 as 



$3 Gj,M x e + x B ') He"e9 a ,9"^e sn ( M e _ u e") 
where the functions Gf are defined by 

He"ee , M x 6 ± x 0") 



(40) 



Gi 



GZ 

Z7T- 



IL» e e Q M x e± x e») ' 
The last factor in (I4UI) can be rewritten using the relation sn(ue — ug> 



s y sm 



sec 



e. g. formula (4.5) in [T3]. The resulting products of sine functions can be calculated 
explicitly: for 9 G 6 a , one has 

9-9" 



1 n 



sm 



-1 Vs 



sm7rz/, 
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2*-i n sin i^! = (_i)^-^-iAT. 

e"ee a ,e"jte 

The remaining products of sinh's can be combined into the function 
FT sinh 7fl+7e " 



n^ e0a sinh^' 

± 



There exists a simple combination of G e and independent of 9, namely 

rt+rt— r^2 

Orn Orn It 



(41) 



with G n = G^. The identity (I4ip can be proven using the standard addition formulas 
for theta functions, their relation to the elliptic functions and the formula (4.8) in [15] . 
relating elliptic and trigonometric parametrization. 

Summarizing the above transformations and using ( T4TT) . we can rewrite the formula 
fl39|) in the following equivalent form: 



1 V AVsinh 7e sinh 7e , tf 3 # 3 (^a' - X a )tf 4 (^ + *v 



(42) 



where the function rjg is given by 

^ = 21n^-ln^. (43) 

Analogous manipulations with (|371)-(l38l) lead to the representations 

detD = (^n MXo -XJ j-j e _ ve/ , j-j ^ (M) 

3 ee6» Q 0e6> , 



D _ x = s y sin Trzy e fa"^')/ 2 fl 2 fl 4 i? 3 (*q' - X Q + x e - x ev 
" AVsinh 7e sinh 7e , tf 3 - X a )A^e - V ' 

By (123]) . the formula ( 144]) gives the vacuum expectation value of twist field. Non-zero 
crossed and non-crossed two-particle form factors are determined by (126]) . (]2Tj) . (145]) and 
(125]). (130]). (P]). respectively. 

Multiparticle form factors f]T6]) can be obtained from ([TS]) — (122]) . They do not vanish 
only if the charge preserving condition m — m' = n — n' is satisfied, where m = #(0), 
n = #(0), vnl = #(#') and n' = #(</>')• Under this restriction, one has 

V ^3 / ^ v / ^smh 7e ^± ^iVsinh^ 

x — a — n /N ■ h n /ft • det ^' ^ 

ViVsmh 7 ,, A g X v/iVsinh^ 

' A B 

where the (m + n') x (m + n') matrix 1Z has 2x2 block form, 1Z = I „ ^ I , and is 



explicitly given by 

i-d 2 {X a i - X a + x dj + Xfa) 
# 3 (X a / - X Q )tf 4 (x + X0 
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$ 3 (X a ,-X a + x 6j -xg>) 
B jk = -^-j^ y n Q / r> 3 = l,---,m, k = l,...,m, 

$z{X a , - X a + x^ h - xq) 

•&z{x a i - x^iix^ - 

i r d 2 {X a - X a , + x 6 ' k + xp.) 
d z (X a - X a ,)$ 4 (x e > + Xtf) 

K J 

The determinant of TZ in (1461) can be evaluated in a closed form using Frobenius 
identity. Given 2L indeterminates z 1 , . . . , zl and z[, . . . , z' L , it expresses the determinant 
of the elliptic Cauchy matrix Q with elements 



T> jk = --rv^? v- n.o ~ — ~ — — » i = . ..,n, fc = l,...,m. 



in the product form 



^(E^-E^+r) n-< fc ^(%-^)^( 



det^ = - ~ J ^±±2^_i Z ' k Z J± . ( 4 J 

Setting L = m + n', r = X a > — X a + | — ^ and 



7TT 

^ + — for j = l,...,m, 



rj — m 



for j = m + 1 , . . . , m + n', 



-x^ for j = 

4 = < vrr 

x 6>' + for j = n + 1 , . . . , n + m , 



one can check that Q coincides with TZ up to diagonal matrix factors. It then follows 
from (iHD that 



•& p (x a > - X a + x e ~ E X6'+J2 x <fi- E x 4>') 

detn = ( _i )m '(n' + i) r ("') 2 ^ ege - e ' ee ' * e * ^ x 

n A(x 6j -xg k ) n ^i(^-^h) n M^-m) n tfiOv fe -^') 

l<j<k<m l<j<k<n l<j<k<m' l<j<k<ri 

m m' n' n 

n n M x oj - x e' k ) n n - ho 

j=l k=l j=l k=l 

m n' n m' 

n n + x <t>' k ) n n + ^) 

J = l k=l j = lk = l 

x- — — ; , 49 

m n n' m' y ' 

n n M x e 3 + x <t> h ) n n + a^j 

j'=l fe=l j=l k=l 

where p = 2 for m - m' e 2Z + 1 and p = 3 for m — m' G 2Z. Together with (T4B]) . this 
gives our main result — a completely explicit factorized formula for any multiparticle 
matrix element of twist field. 
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3.4- Thermodynamic limit 

Form factors found above simplify in the thermodynamic limit N — > oo. To explain 
these simplifications, consider e. g. the expression ( 1441) for the vacuum expectation 
value ( 125|) . We need to evaluate the asymptotics of (i) X a i — X a and (ii) i]g. Both 
quantities can be written in the same form 

e'ee a , e>ee a 

where the function fgr is defined for 6' e [0, 2ir]. In the case (i) one has fgt = xg>, while 
in the case (ii) /#/ is given by 

^' =21n T7 — T r~ lnsmh o • 51 

V4{Xg + Xg>) 2 

The most important distinction between the two situations is that in the first case one 
has f 2n — fo = 7T, whereas (loTj) extends to a continuous 2-7r-periodic function on the real 
line. 

As N — > oo, one has the estimate 

f f 2?rzy f> 

AT-1 2tt 

with & = 0, . . . , N — 1. The sum ^ transforms into the integral f f' g d6, which 

fc=0 N o 

yields 

lim (X a i — X a ) = 7Yb>, lim rjg = 0. (52) 
In the same way one shows that 

a(E»-E*)=°- ( 53 ) 

0e0 a , eeo a 

Rigorous mathematical proofs of fl52l) - fl53l) and explicit error estimates can be obtained 
using a kind of Sommerfeld- Watson transform of the sums ( 1501) . 

It then follows from ([25]) . (|44j) that the thermodynamic limit of the vacuum 
expectation value of twist field is given by the simple expression 

(a) =' lim a {vac\O a>a ,{0*)\vac) a , = ^PT, (54) 

where the arguments of theta functions and the half-period ratio are now indicated 
explicitly for further convenience. As one could expect on general grounds, the r.h.s. of 
(I54p depends on a, a' only via the difference v = a' — a. Note that this relation has the 
same form as the first formula on p. 187 of [25] . 

Particle quasimomenta in the thermodynamic limit uniformly fill the interval [0, 2tt]. 
Summation over each of them in the form factor expansions of correlation functions 
transforms into integration: J2e ~~ 2^ Jo * ^9- Although this naive procedure is 
plagued by the appearance of annihilation poles in the crossed form factors, it suggests 
to consider 

F„(0;0|0';<//) = f lim N ^¥^ jr { 0] 0|0' ; 0'). 

N—>oo 
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The asymptotics (j52p implies that the limiting two-particle form factors are given by 



T v {e-\e'-) = T u {-e\-e>) = 

SySmiiveW- - 2 ^' 2 ^ 2 (0|r)^ 4 (0|r) T&^xg - x e > + nu\r) 
Vsinh 7^ sinh 7^ $ 3 (0\t)$ 3 (ttv\t) i?i (x e - x e > | r) 

F,(^^|;) = -F_(;|^^) = 

is^sinvrz/e^ 9 '- )/ 2 ^ 2 (0|r)^ 4 (0|r) # 2 (x0 + x fl / + irv\r) 



(55) 



(56) 



Vsinh 70 sinh 79/ ^3(0|r)t? 3 (7rz/|r) # 4 (x0 + xg>\r) 
The limit of multiparticle form factors can be found in a similar way: it suffices to 
remove the functions rjg from ( 146]) and to replace X a > — X a in (1491) by 7tza 
In the case f = ±~, the formulas ( I5"4"j) - (156"j) reduce to simpler expressions 

{0 ± ,) = {l-k 2 )\ (57) 

ie i(e'-S)/2 ginh Tfl+Tfii 

F ± i(0; 6K; ) = , e _ 2 , , (58) 

2 A/smh 7 e sinh 7^/ sin 

-is*s e i{e '- 9)/2 sin^ 
F . ! (0; ) = / W 2 . (59) 

2 V sinh 76) sinh 79/ sinh 79 

They can be proved using the relations #1(2+!) = $2(2), $3(2+2) — $4(2)1 rewriting the 
theta function ratios in terms of elliptic functions and finally passing to the trigonometric 
parametrization. 

3.5. Scaling limit 

The gap in the energy spectrum closes as 70 — > 0. Since 70 = 2(/C y — /C*), this 
corresponds to k — > 1 (or, equivalently, r — >■ zO or q — >■ 1) in terms of elliptic parameters. 
Also note that 

-t- x = — = — In— — + o(l) as&^l. 60 
K ix 16 

In order to obtain the asymptotics of the vacuum expectation value (1541) in the vicinity 
of the critical point, recall the modular transformation 

.2 



■& 3 {z\t) = {-ir)-y 2 e-^^[--\ - -J, (61) 
and the product formula 

3 ( Z \ T ) = Y[{1 - q 23 )(l + q 23 ~ x e 2iz ){\ + q 2 ^ l e~ 2iz ). (62) 
i=i 

Rewrite the theta functions in (154")) using (l6Tj) . The nome of the transformed functions 
vanishes as k — > 1. The representation (16"2"|) therefore implies that 



/ 7rz/, 1\ [1 for I r/| < \ 

limtf 3 --) = <^ 1 

r^io V r r/ 2 for 1/ = ±± 
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This in turn can be used together with (1601) to show that, as k — > 1, 

(a)= (2-- (l - fc T' [1 + o (1 )] fa |„|<i, m 
{ (l-/c 2 ) 3 [l + o(l)] fori/ = ±i. 

The second approximation is in fact exact, cf. (1571) . As the r.h.s. of (|54|) is periodic 
in i/ with period 1, the critical asymptotics of the vacuum expectation value for any v 
can be deduced from fl63|) . Note that the scaling dimension of the twist field for \v\ < | 
has the expected value v 2 . 

In the vicinity of the critical point, our initial lattice model becomes equivalent 
to a field theory of free massive Dirac fermions. Correlation functions of twist fields 
in this theory are determined by momenta at the scale of inverse correlation length 
(for convenience, we change the domain of definition of lattice momenta from [0, 2ir] to 

[-7T,7T]). 

1 - k 2 

More formally, denote e = , set 9 = £sinh£ and let k — > 1. The dispersion 

2c y 

relation Q then implies that 70 — > es y cosh£. One also has 

Xg - Xqi _^ i(£ - £') 



T 2 

Normalized scaled two-particle form factors of twist fields are determined by 
F„(£;|£';) = limeF I/ (esmh£;|£sinh£ , ;) > 
F„(£;£'|;) = hm £ F„( £ smh£;-£sinh £'];). 

k— >1 

where the variables £, £' have the meaning of particle rapidities. 

To compute the corresponding limits, one can adopt the same approach as in the 
above asymptotic analysis of the VEV. Transform the theta functions in (j55p - ()56p using 
Jacobi's imaginary transformations, rewrite the result using product formulas, and then 
let r — > iO. For \v\ < |, one finds 

sin ttu pKf'-f-*"') 

F„(f ; |f; ) = — - — , 

Vcosh | cosh I 7 sinh 

F (C'C'I - ) siniru ie v ^~® 

Vcosh | cosh £' cosh 

This reproduces two-particle (and hence all) form factors of the continuum twist fields 
in the massive Dirac theory (TJ [23l EHJ El], which correspond to the exponential fields 
in the sine-Gordon model at the free-fermion point. 

Likewise, it can be deduced from (I58l) - (l59l) that for v = ±~ one has 

— i coth ^—^ 



F ±i (£;!£';) = ^[F|_ (e;ie';) + F_i +0 (e;ie'; 



2 



Vcosh^ cosh 



tanh ^rr^ 



2 



y/ cosh I cosh 
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The fact that Fi_ 7^ 1F_i+o can be understood as follows. Monodromy conditions for 
fermion fields in the continuum give a system of integral equations for the two-particle 
form factors pQ. This system has a unique admissible solution for \v\ < |, and two 
solutions for v = ±|. Each of the two solutions gives rise to a twist field operator. Any 
linear combination of these operators leads to the required fermion branching. 



4. Concluding remarks 

Finite-lattice form factors in the conventional free-fermion models (triangular Ising 
lattice, XY quantum spin chain, BBS2 model) can be obtained [121 E3] from those 
of the Ising spin on the square lattice [31 HI El QUI C3]. The fields considered in this 
paper are more general; in particular, their scaling dimension continuously depends on 
a real parameter v. We compute their form factors explicitly in terms of the Jacobi 
theta functions and show that in the scaling limit they reduce to form factors of the 
exponential fields of the sine-Gordon model at the free-fermion point. 

Determinantal form of form factors is also encountered in some of the interacting 
integrable models, see e. g. Slavnov's formula for scalar products of Bethe states in the 
spin-i XXZ chain [T71 [32]. An intriguing question is therefore if it is possible to go 
beyond the free-fermion point and calculate form factors of twist fields in the integrable 
lattice regularization of the massive Thirring model [22] (related to the eight-vertex 
statistical model). We believe that the present work makes a step in this direction. 

Another challenge is to complete form factor derivation for Z^r-symmetric 
superintegrable chiral Potts quantum chain. It was recently shown that they have Ising 
form up to unknown scalar factors labeled by the pairs of the so-called Onsager sectors 
[T5] . Putative free-fermion part of the integrable structure of this model is yet to be 
elucidated. 
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